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Bessel Functions of Real Argument and Integer Order 
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A computer program is described for calculating Bessel Functions J n (x) and l H {x), for x real, and 
n a nonnegative integer. The method used is that of backward recursion, with strict control of error. 
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1 . Method 

Given a real number .% and a positive integer NB, BESLRI calculates either 

/„(*), n = 0,l,. . .,NB-1 

or 

,/„(%), n = 0,1,. . .,NB-1 

using double-precision arithmetic. The method, which is described in [l], 2 is based on algorithms of 
Olver [2] and Miller [3], applied to the difference equation 

In 
Jn-\ =— Jn~ SIGN • y n +i 

where SIGN is ■+ 1 for J 's, — 1 for/ 's. 

The program sets MAGX= [ | x | ], the integer part of | x |, Pmagx = 0, Pmagx +i = I> and 
then successively calculates 

p w+ i = f^TPn "SIGN • pn-i zi = MAGX+ 1, MAGX + 2, . . .. (2) 

The sequence is strictly increasing. The program takes N to be the least n such thatp /t exceeds a 
number TEST defined in sections 2 and 3. It then sets y[ N) = 0\ y { f_ x =Hp.x, and recurs backward 
using (1). The computed sequence y { o\yi N \ ... is the recessive solution of (1) which satisfies the 
boundary condition ^magx = 1. From this solution, the /'s and 7 's are found by normalizing 

Jn(x)=yi^lfi n = 0, 1, . . .,NB-1 
In(x)=yWlfJL n = 0, 1, . . .,NB-1 
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M = y ( V v, + 2 2 2l y^ fori's 

k= 1 



li = (y,V V) + 2 [ ff ] y$)/cosh x for /'s. (3) 



2. Error Bounds for the y n 

For /i > MAGX, the truncation error in y { n N) is 



^ PrPr+l ' 



see [2], eqs (5.01) and (5.02). This error is bounded by using Lemma 2 of [1], which states that for 
s^r 9 p s = p s+1 lp s ^ min (Pr+i/Pr, (r+l)/|*| + V(r+ 1) 2 /|%| 2 - 1). 
The program sets 



TEST, ^ V2-10Nsi« PaP ^ +1 , (4 ) 

where L = max (MAGX -h 1, NB — 1), and NSIG is the maximum number of significant decimal 
digits in a double-precision variable on the computer being used. Then N' is the least n such that 
p n > TESTi. For/ 's, N is the least n ^/V' such that 



p n > TEST = J-^t- TEST,. (5) 

V pf - 1 

In consequence of (4) and (5), the relative truncation error | T { fj ) / y„ | is less than - • 10~ NSI(; for all n 

in the range MAGX < n ^ L; see [1, sec. 5]. For /'s, /V is taken to be N'. Here a glance at (1) shows 
Pn I Pn- -i > 2 for all ai > MAGX, so the proof given in [1, sec. 5] for J's may also be used for /'s, 
with p.v replaced by 2 throughout. 

In eq (1) for /'s, — y ( ^ } and —SIGN • yW have the same sign since if x > 0, then yW > Ofor 

all m, and if x < 0, then consecutive y^'s have opposite signs. Thus relative errors < e in y^/ and 
y ( JP will produce a relative error < e in y\ t N _\. Therefore, the relative truncation error in all I n (x), 

^ x < NB, is bounded by ^ • 10 NSI(: . 

For J n (x) , n ^ MAGX, the relative truncation error cannot be bounded, since as n = MAGX, 
MAGX — 1, . . ., 1 in (1), the values J n (x) oscillate, and precision is lost owing to cancellation. In 
this range, J n (x) is accurate to about D decimal places, where D is the number of decimals in 
J MM.x+i 0*0 which corresponds to NSIG significant figures in the same quantity [1, sec. 5j. 

3. Normalization and Error Bounds for fx 

The eqs (3) were chosen to keep cancellation under control. First, 

\Mx) |*1, \Jn(x) l<^=. forn^l (6) 
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[4, 2.5], so each term of the sum (3) for 7's is less than v2 times the whole sum. Since (6) is a 
rather weak bound, cancellation is even less than this would indicate. 

For /'s, cancellation is avoided altogether, since all terms in the sum (3) for/'s have the same 
sign, as shown in section 2. 

Besides bounding the truncation error of the algorithm, the program provides an estimated 
bound for the truncation error of the normalization sum, defined by eqs (3). For/'s, this error is 



S<")=r - r (">+2 £ (y2 k -y^ 



For 2A:^ MAGX, a bound for the error term y^-y-ik is unavailable. For MAGX < 2k < /V, 
I y\V~ y*k I < P2kP.\l(p\ — 1) [1, sec. 5]. To avoid storing all the p 2 k, the program allows only 

N 

for terms for which 2k ^ N. Here }il ) = 0, and y r 2k = Pik V l/(p r Pr + i ). Therefore 



yijp-yu 



Pik + i 



P2k Pik />2frfl 

P2k+2 P2kr> P'2k+3 



P2k+1 



1 1 

1 + + - -+. . .\ 

plk pik 



p%> 



(p*., - i )p 2 fe+i 



compare section 2 above. Now let 



The 



R ( " ] = 2 \y» 



*- [-i-] 



/?(" 



P 2 v 



P 2 v 



1 ■" Pik f i 

r -v + i t 



(p-y, - 1)P.V+1 



P2, P 4 V , 



P 3 v 



(7) 



(p 2 v ,-l) 2 p.v 
ForV's the program sets TEST i 3= 2 • 10 NSI(; . The normalization factor p is 1IJmm;x(x) [1, sec. 5], 



R(») 



Jmai;x(-*') 



R(") 



./ma<;x(*) 1 P 3 v . 1 
(p.v'-l)* 2 



10" 



This is a rather weak upper bound, and the error 

For / 's, Pnlpn - 1 > 2 so the analog of (7) is 



S (.v) 



p 



turns out to be less than-- 10 Ns "' 



fl(.v) <s 



8 



9p,v cosh x' 
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this is derived by substituting 2 for p.v in (7), and introducing the factor cosh x; see (3). 

2- io -nsig 
By setting TESTi ^ — ZtFTJZi \/t\rv\ > tne P r °g ram insures that 



1 ^ exp(0.461 • MAGX)' 

( v) < 4 exp(0.461 - MAGX) 



9-10- NS,(: -cosh 



Here, jul is -: (x), so the relative error | 7? ( v) /m I is bounded by 

JMACX 



#(Y) 



/* 



4 /magx(*) exp(0.461 MAGX) 
9 cosh x - i0- NSI(i 



(8) 



Now Kapteyn's inequality [4, 8.7] states that 

z exp Vl— z l 



\Jn(nz) 



1+ Vl-2 2 



Let z 



Then 



M . r Y , so for large *, z — 1, approximately. 



|/magx(*)I = |/magx(**)I = 1 7 magx (MAGX • iz)\ 



z exp Vl +_z 2 
1+ Vl+z 2 



For large %, the approximate bound for /magxM is 



exo V? \ MAGX 

y^| ) < exp (0.533 ■ MAGX). 



Substituting this in (8), we obtain 



RW 
^ 



< 



4 • exp (0.533 MAGX) ; exp (0.461 ; MAG X) 
9 cosh x • 10" NSIG 



< - • 10- NSIG < - • io- nsig . 



Besides this approximate upper bound, the strict bounds 






< 7.6 • 10 



-NSIG 



and 






< 0.565 • 10 



x < 1 



may be derived by using Kapteyn's inequality to obtain a strict bound for |/magx+i W 
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4. Appendix: Algorithm BESLRI 



ELT BESLRI* 1*730222* 5*730 



000001 
000002 
000003 
OO0O04 
000005 
000006 
000007 

ooooob 

000009 
000010 
000011 
000012 
000013 
000014 
000015 
000016 
000017 
000018 
000019 
000020 
000021 
000022 
000023 
U00024 
000025 
00002b 
000027 
000028 
00U029 
000030 
000U31 
000032 
000033 
000034 
00003b 
00U036 
000037 
000038 
000039 
000040 
000041 
000042 
000043 
000044 
000045 
00004b 
000047 
C00048 
000049 
000050 
000051 
000052 
000053 
000054 



C 
C 
C 
C 
C 
C 

c 
c 
c 
c 
c 
c 
c 

C B 



SUBROUTINE BESLRI ( X . NB » I ZE* B * MC AI.C ) 
THIS ROUTINE CALCULATES BESSEL FUNCTIONS I AND J OF REAL 
ARGUMENT AND INTEGER ORDER. 



EXPLANATION OF VARIABLES IN THE CALLING SEQUENCE 



NB 
IZE 



DOUBLE PRECISION HEAL ARGUMENT FOR WHICH I*S OR J*5 
ARE TO BE CALCULATED. IF I*S ARE TO RE CALCULATED* 
ABS(X) MUST BE LESS THAN EXPARG (WHICH SEE BELOW). 
INTEGER TYPE. 1 + HIGHEST ORDER TO BE CALCULATED. 
IT MUST BE POSITIVE. 

INTEGER TYPE. ZERO IF J*S ARE TO BE' CALCULATED * 1 
IF I*S ARE TO BE CALCULATED. 

DOUBLE PRECISION VECTOR OF LENGTH NR» NEED NOT BE 
INITIALIZED BY USER. IF THE ROUTINE TERMINATES 
NORMALLY (MCALC=NB)r IT . RETURNS J (OR D-SUR-ZERO 
THROUGH J(OR I ) -SUB-NB-MINUS-ONE OF X IN THIS 
VECTOR. 
NCALC INTEGER TYPE* NEED NOT BE INITIALIZED BY USER. 

BEFORE USING THE RESULTS* THE USER SHOULD CHECK THAT 
NCALC=NB» I.E. ALL ORDERS HAVE BEEN CALCULATED TO 
THE DESIRED ACCURACY. SEE ERROR RETURNS BELOw. 



EXPLANATION OF MACHINE-DEPENDENT CONSTANTS 

G DECIMAL SIGNIFICANCE DESIRED. SHOULD RE SET TO 

IFIX(ALOGin(2)*NBIT+l ) r WHERE NBlT IS THE NUMBER OF 
BITS IN THE MANTISSA OF A DOUBLE PRECISION VARIABLE. 
SETTING NbIG LOWER WILL RESULT IN DECREASED ACCURACY 
WHILE SETTING NSIG HIGHER WILL INCREASE CPU TIME 
WITHOUT INCREASING ACCURACY. THE TRUNCATION ERROR 
IS LIMITED TO T=.5*10**-NSIG FOR J*S OF ORDER LESS 
THAN ARGUMENT* AND TO A RELATIVE ERROR OF T FOR 
I*S AN[) THE OTHER J.*S. 
NTE K J LARGEST INTEGER K SUCH THAT 10 + + K IS MACHINE- 
REPRESENTABLE IN DOUBLE PRECISION. 



C LARGEX UPPER LIMIT ON THE MAGNITUDE OF X. BEAR IN MIND 



THAT IF ABS(X)=Nr THEN AT LEAST N ITERATIONS OF THE 
BACKWARD RECURSION WILL BE EXECUTED. 
EXPARG LARGEST DOUBLE PRECISION ARGUMENT THAT THE LIBRARY 
DEXP ROUTINE CAN HANDLE. 



ERROR RETURNS 

LET G DENOTE EITHER I OR J. 

IN CASE OF AN ERROR* NCALC. NE.NR* AND NOT ALL G*b 
ARE CALCULATED TO THE DESIRED ACCURACY. 

IF NCALC. LT.Oi AN ARGUMENT IS OUT OF RANGE. NB.LL.O 
OR IZE IS NEITHER NOR 1 OR IZE=1 AND ABS ( X ). GF .EXPARG. 
IN THIS CASE' THE B-VECTOR IS NOT CALCULATED* AND NCALC 
IS SET TO MINO(NB*0)-1 SO NCALC. NE.NB. 



1 
00000200 
00000300 
00 0004 GO 
00000500 
00000600 
00000700 
00 H 
00000900 
0000 1000 
00001100 
00001200 
00U01300 
00001400 
00001500 
00001600 
00001700 
00001800 
00001900 
00002000 
00002100 
00002200 
00002300 
00002400 
25 
00002b()0 
00002700 
00002800 
00002900 
00003000 
00003100 
00003200 
00003300 
00003400 
00003500 
00003600 
00003700 
00003800 
OOU03900 
00 04 00 
00004100 
00004200 
00004300 
u 4 4 
00004500 
000046 
00004700 
001104600 
00004900 
00U05000 
00005100 
00005200 
00005300 
00005400 
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1)00055 C NB.GT.NCALC.GT.O WILL OCCUR IF NB.GT.MAGX AND ABS(G- 00005500 

000056 C SU3-NB-OF-X/G-SUB-MAGX-OF-X).LT.10.**(NTE"N/2) » I. E. NB 00005600 

U00057 C IS VUCH GREATER THAN M A GX . IN THIS CASE> R(N) IS CALCU- 00005700 

UOU059 C LATED TO THE DESIRED ACCURACY FOR N.LE.NCALC* BUT FOR 00005800 

U00059 C NCALC.LT.N.LE.Nn, PRECISION IS LOST. IF N.GT.NCALC AND 00005900 

U00060 C ABS<B(NCALC>/B(N)).EQ.10**-K» THEN OMLY THE FIRST NSIG-K 00006000 

U0J061 C SIGNIFICANT FIGURES OF B(N) MAY RE TRUSTED. IF THE USER 00006100 

000062 C WISHES TO CALCULATE B(N) TO HIGHER ACCURACYr HE SHOULD USE 00006200 

000063 C AN ASYMPTOTIC FORMULA FOR LARGE ORDER. 00006300 
UO.JUbU C 00006400 
000065 DOUBLE PRECISION 00006500 
U00u6o 1 X#rf»P»TEST»TEMPA#TEV!P3»TEMPC»EXPARG»SlGN»SUM#TOVER» 00006600 
U00l)b7 2 PLAST»POLD»PSAVE#PSAVEL 00006700 
000060 DIMENSION B(MB) 0OO06BO0 
O0liib ; J DATA MS1G»NTFN»LAR6EX»EXPARG/19»307»100000»7.D2/ 00006900 
O0D070 TEMP£=DABS(X) OOU07000 
U0OU71 MAGX=IFIV(SNGL(TEMPA)) 00007100 
000 72 IFd>JB.GT.0.AMO.MAGX.LE.LARGEX.AND. ( I7E.EQ . 0.0R. 0000 72 00 
000073 1 ( IZE.EQ. I . ANH.TEMPA .LE.EXPARG) ) ) GO TO 1 00007300 
00007-4 C E«RO« RETURN -- X . NB * OR IZE IS OUT OF RANGE 000071*00 

000075 NCALC=MlNO(NB»0)-l • 00007500 

000076 RETURN 00007600 

000077 1 SIGN=DaLE(FLOAT(l-2*IZE> J 00007700 
GOU07B NCALC=NB 00007800 
001)079 C USE 2-TERM ASCENDING SERIFS FOR SMALL X 00007900 

000080 IF(TEMPA**4.LT. .1D0**NSIG) GO TO 30 00008000 

000081 C INITIALIZE THE CALCULATION OF P*S 00008100 

000082 MBV»X=N8-MAGX 00008200 

000083 N=MAGX+1 00008300 
00008'* PLASTrl.DO 00008400 

000085 P=DJLE(FLOAT(2*N) J/TEMPA 00008500 

000086 C CALCULATE GENERAL SIGNIFICANCE TEST 00008600 

000087 TEST=2.D0*1.D1**NSIG 00008700 

000088 • IF( IZE.Eu.l .AND.2*MAGX.GT.5*NSIG) TEST=DSQRT ( TEST*P ) 00008800 
00089 IF( IZE.EQ. LAND. 2 *MAGX.LE.5*NSIG)TEST=TEST/1.585**MAGX 00 890 

000090 'A-0 0OU090O0 

000091 IF(NRMX.Lf .3) GO TO 4 00009100 
OOU092 C CALCULATE P*S UNTIL M=NB-l. CHECK FOR POSSIBLE OVERFLOW. 00009200 

000093 TOVEH = 1.01**(NTEN-,MSIG) 00009300 

000094 NSTART=MAGX+? - 00009400 
00U095 IMEND=N8-1 00009500 
00u09o DO 3 N=NSTARTrNEND 00009600 

000097 P0L0=PLAST 00009700 

000098 PLAST=P , 00009800 
00u09<i P=DRLE (FLOAT (2*N) )*PLAST/TEMPA-SIGN*POLO 00009900 

000100 IF(P-TOVER) 3»3»5 00010000 

000101 3 CONTINUE 00010100 

000102 N=NEND 00010200 

000103 C CALCULATE SPECIAL SIGNIFICANCE TEST FOR NRMX.GT.2. 00010300 
00104 TEST=DMAX1(TEST»DSQRT(PLAST*1.D1**NSIG)*DSQRT(2.D0*P) ) 00010400 

000105 C CALCULATE P*S UNTIL SIGNIFICANCE TEST PASSES 00010500 

000106 4 N=N+1 00010600 

000107 POLD=PLAST 00010700 
OOOlOy PLA5T=P 00010800 

000109 P=DBLE (FLOAT ( 2*N) )*PLAST/TEMPA-SIGN*POLD 00010900 

000110 IF(P.LT.TEST) GO TO 4 00011000 

000111 IF(IZE.EO.l.OR.M.EQ.l) GO TO 12 00011100 

000112 C FOR U*S» A STRONG VARIANT OF THE TEST IS NECESSARY. 00011200 
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000113 C CALCULATE IT, AND CALCULATE P*S UNTIL THIS TEST IS PASSEU. 00011300 

000114 M=l 00011400 
00011b TEMPB=P/PLA5T OoollbOO 
000116 TE^PCrOBLE ( FLOAT (N+DJ/TEMPA 00011600 
0(1117 IF( TEMpB»l.Dn/TE | WIPB.GT.2.D0*TEMPC)TEMPR=TEMPC+DSQRT 00U11700 

000118 KTEMPC**2-l.nO) 00011800 

000119 TEST = TEST/DSQRT(TEMPB-1.00/TEMPD) 00011900 

000120 IF(P-TEST) 4»12»12 00012000 

000121 C TO AVOID OVERFLOW, DIVIDE P*S BY TOVER. CALCULATE P*b 00012100 

000122 C UNTIL ABS(P).GT.l. , 00012200 

000123 5 T0VER=1.D1**NTEN 00012300 

000124 PrP/TOVER 00012400 
00012b PLAST = PLA5T/T0VER 00012500 

000126 PSAVE^P 00012600 

000127 PSAVEL=PLAST 00012700 

000128 NSTART=N+1 00012800 

000129 6 N=N+i ' 00012900 

000130 POLD=PLAST 00013000 
0UU131 PLAST=P 00U13100 

000132 P=DULE(FLOAT(2*N))*PLAST/TEMPA-SIGN*POLD 00U13200 

000133 IF(P.LE.l.DO) GO TO ft 00013300 

000134 TEV,PB = DBLE(FL0AT(2*N) )/TEMPA 00013400 
00013b IF(lZF.EQ.l) GO TO 8 00013500 
00013b TEVPC=.5D0*TEMPB 00013600 

000137 TEV|PR=PLAST/POLO 00013700 

000138 IF(TEMPn+l.Dn/TEMPB.GT.2.D0*TEMPC)TEMPB=TEMPC+DSQRT ' 00013H00 

000139 L(TEMPC**2-l.O0) 00013900 

000140 C CALCULATF BACKWARD TEST, AND FIND NCALCf THE HIGHEST N 00014000 

000141 C SUCH THAT THE TEST IS PASSED. 00014100 

000142 8 TEST=.bD0*POLD*PLAST*(l.D0-l.D0/TEyPB**2)/l.Dl**NSlG 00014200 

000143 P=PLAST*TOVER 00014300 

000144 N = N-1 00014400 
00014b NENl)=MIN()(Nli»N) 000l4b00 

000146 DO 9 NCALC=NSTART»NEND 00014800 

000147 POLD=PSAVEL 00014700 

000148 PSAVELrPSAVE 00O14MOO 

000149 PSAVE = n»LE(FL0AT(2*N) ) *PSAVEL/TEMPA-SIr,N*POLD 00014QQO 

000150 IF(PSAVE*PSAVEL-TEST) 9»9»10 000151)00 
OOOlbl 9 CONTINUE OOOlbluO 
0001b2 NCALC-NEND+1 . 000lb200 

000153 10 NCALC=NCALC-1 00()lb300 

000154 C THE SUM 3 ( 1) +2B(3)+2B<5) • • « IS USED TO NORMALIZE. M» ThL 00015400 
OOOlbb C COEFFICIENT OF B(N), IS INITIALIZED TO 2 OR 0. 00015500 

000156 12 N=N+1 OOOlbf.OO 

000157 M=2*N-4*(N/2) 00015700 

000158 C INITIALIZE THE BACKWARD RECURSION AND THE NO&MALIZA T ION 00015800 

000159 C SUM 01)015900 
000180 TEMPi3=0.D0 001)16000 

000161 TEMPA=1.D0/P OOOlblOO 

000162 SUM=DHLE(FLOAT(M) )*TEMPA 00016200 

000163 NEND=N-Nf3 00016300 

000164 IF(NFNJ) 17,15,13 001)16400 
00016b C RECUR BACKWARD VIA DIFFERENCE EQUATION, CALCULATING (BUT 00016500 

000166 C NOT STORING) B(N), UNTIL M=NB. 000lt>600 

000167 13 DO 14 L=lfNEND 00016700 

000168 N=N-1 00016800 

000169 TE^PC^TEMPB 00016900 

000170 TE'1PB = TEMPA 00017000 
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000171 
0O0172 
000173 
000 17U 

0017b 
00 17b 
000177 

o o u 1 7 a 

17 9 
000180 

o o o i b i 

000132 

1 8 3 

Oil 184 
13b 
UOOlfao 

00 187 

uooia* 

00 189 
00D190 
U 1 9 1 
00192 
U00193 
00 194 
00 19b 
00019b 
000 197 
00019H 
000199 
0l)u200 
000201 
00202 
0u20 3 
00 20'+ 
0020b 
00206 
2 7 
000 20 8 
2 U 9 
210 
211 
UUJ212 
213 
214 
00 21b 
00216 
217 
00 218 
00021^ 
000220 
000221 
0222 
3. 



TEMPA=(DBLF( FLOAT (?*N) ) *TEMPB ) /X-STGN*TEVPC 
M=2-M 

14 SLn=S'JM+DBLE< FLOAT (M) ) +TEMPA 
C STORE B<NR) 

15 B(N)=TEMPA 
IF(NM.GT.l) GO TO 16 

c nb=i. since 2*TEMPA was added to the sum, tfmpa must be 
c subtracts 

sum=sum-tempa 

GO TO 25 
C CALCULATE AND STORE B(NH-l) 
lb M=-M-l 

B(N) -(D n ,LE(FLOAT(2*N) ) *TEMPA ) /X-SIGN' + TF^PB 

IF(N.EQ.l) GO TO 22 

M=2-M 

SUM=SUM+DBLE (FLOAT (M) ) *B(N) 

60 TO 19 
C N.LT.NB* SO ST0R r B(N) AND SFT HIGHER ORDERS TO ZERO 

17 t3(N)=TE^PA 
NEN,J = -NEND 

00 IP, L-l »NEND 

18 B(M+L>=0.Dn 

19 MEMO=N-2 
TF(UENO.FO.O) GO TO 21 

C CALCULATE VIA DIFFERENCE EQUATION AND STORE B(N)» 
C UNTIL N=? 

DO 20 L=l »NE'ID 

N=N-1 

B('i) = (DBLE (FLOAT (2*N) )*B(M+1) ) /V-SIGN*" ( N+2 ) 

M=?-M 

20 SUV=SUM+DBLE< FLOAT (M) )*B(N) 
C CALCULATE 8 ( 1 ) 

21 3< 1)=2.00*B(2)/X-SIGM*B(3) 

22 SUM=SUM+B<1) 

C NORyALIZF— IF IZE=li DIVIDE SUM BY COSH(X). DIVIDE ALL 
C (j(N) BY SUM. 

23 IF( IZE.EO.O) GO TO 2b 
TEMPA=DEXP(DABS(X) ) 
SUM=?.DO*SUM/(TEMPA-»-l.D0/TEMPA) 

2b DO 26 N=l »NB 
2b 9(U)=B(N)/SUM 
RETURN 
C TWO-TERM ASCENDING SERIES FOR SMALL X 

30 TEMPArl.nO 
TEMPR=-.25D0*X*X*SIGN 
tt< 1 )=1.00+TE M PR 
IPCJiJ.EQ.l ) GO TO 32 
DO 31 N=2»NB 

TE v f , A-TEVPA*X/DBLE(P-L0AT(2*N-2) ) 

31 B(N)rTEMPA*(l.D0+TEMPB/D3LE(FLOAT(M) ) ) 

32 RETURN / 
END 



00017100 
00017200 
00017300 
00017400 
00017500 
00017600 
00O17700 
00017H00 
00017900 

oooiaooo 

01810 
00 010200 
00018300 
00 0184 00 
00018500 
00010600 
00018700 

oouioaoo 

000l8 r >00 
0OU19O00 
00019100 
00 192 00 
00019300 
00 01 90 00 
00019500 
00019600 
00019700 
00019800 
00019900 
00020000 
00020100 
00020200 
00020300 
00020400 
00020500 
00020600 
00020700 
00020800 
00020900 
00021000 
00021100 
00021200 
00021300 
00021400 
00021501 
00021600 
00021700 
00021800 
00021900 
00022000 
00022100 
00022200 



END CUR LCC U02-0038 L8 
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